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Abstract 

It is shown that the method of lowest order constrained variational (LOCV) 
which is based on the cluster expansion theory is a reliable many-body tech- 
nique to calculate the nuclear matter equation of state. In this respect, the 
state dependent correlation functions and the effective interactions which have 
been produced by the LOCV calculation with the Reid soft core and the A- 
Reid interactions are used to estimate the size of higher order cluster terms 
such as the effect of three-body cluster energy on nuclear matter. We find that 
the three-body cluster energy is less than 1 MeV beyond the nuclear matter 
saturation density and it has weaker density dependence than our previous 
calculation with the state-averaged correlation function and the effective in- 
teraction. Finally we conclude that the LOCV method is good enough to 
calculate other properties of quantal fluids. 
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I. INTRODUCTION 



The method of lowest constrained variational (LOCV) was developed in 1977-1979 [1,2] to 
calculate the bulk properties of homogeneous nuclear fluid such as the saturation energy and 
density, the surface energy and the asymmetrical coefficient in the semi-empirical formula, 
by using the realistic nucleon-nucleon interactions i.e. the Reid [3] and the A-Reid [2] 
potentials. The LOCV method was successfully extended further for finite temperature 
calculation [4] and various thermodynamics properties of nuclear matter were evaluated. 
Recently this approach was reformulated to include more sophisticated interactions such as 
the Vu [5], the AVu and the new Argonne AVis [6] potentials. In these calculations [5,6] 
it was shown that the new AVis potential, like the other N-N interactions, dose over bind 
nuclear matter and a very good agreement was found between the LOCV results and more 
complicated approaches such as the Bruckner-Hartree- Fock (BHF), the Correlated Basis 
Function (CBF), the Bruckner-Bethe (BB) and the Variational Hypernetted Chain (VHC) 
techniques. The results of these calculations are demonstrated in table 1 by presenting 
the saturation energy, the saturation density and the incompressibility of nuclear matter 
for different potentials and different many-body methods from reference [6]. According to 
this table, it is seen that only the results of the AVu and the A-Reid potentials are close 
to the empirical values and on average the inclusion of density dependence three-nucleon 
interaction (TNI) improves the calculations. 

In order to test the convergence of our LOCV results for nuclear matter, we performed 
calculations beyond the lowest order (the two-body cluster term) and the three-body cluster 
energy was calculated with the state-averaged correlation function [10] which in turn was 
generated by using the state-dependent LOCV correlation functions. The smallness of the 
normalization parameter (the convergence parameter) and the three-body cluster energy 
where indicated that at least up to the twice nuclear matter density, our expansion converges 
reasonably and it is a good approximation to stop after the two-body cluster terms [2,5,10]. 
LOCV formalism has several advantages with respect to the other many-body techniques like 
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the Brueckner-Bethe (BB) [7] , the Variational Hypernetted Chain (VHC) [8] , the Correlated 
Basis Function (CBF) [8] and the Brueckner-Hartree-Fock (BHF) [9] which go beyond the 
lowest order: 

1) There is no free parameter in the LOCV method beside the N-N potentials i.e. it is 
fully self-consistent. 

2) It considers constraint in the form normalization condition [11] to keep the higher 
cluster terms as small as possible. As we pointed above, this has been tested by calculating 
the three-body cluster terms with the state-averaged correlation functions. 

3) It assumes a particular form for the long-range part of the correlation functions in 
order to perform an exact functional minimization of the two-body energy with respect to 
the short-range behavior of the correlation functions. 4) It dose functional minimization, 
rather than simply parameterizing the short-range behavior of the correlation functions. So 
in this respect it also saves an enormous computational time. 

The aim of present work is to calculate the three-body cluster terms with the state- 
dependent correlations i.e. without making the state-averaged approximation. In this work 
we use the Reid soft core [3] and the A-Reid [2] potentials, since we are only interested in the 
size of three-body cluster energy. Beside this the results of our previous calculation shows 
that the A-Reid interaction can reasonably reproduce the properties of nuclear matter. On 
the other hand the A state is the most important configuration which modifies the nuclear 
forces and one may consider it as the origin of understanding of the three-body forces [12]. 

So the paper is organized as follows. We begin section II by describing briefly the LOCV 
formalism. Section III-A is devoted to the three-body cluster energy and the definition of the 
state-averaged correlation function and the effective interaction. In section III-B we present 
the explicit equations for the three-body cluster energy in terms of the state- dependent 
correlation functions and the effective interactions. Finally, in the last section we discuss 
our results, the various aspects of LOCV formahsm, the three-body cluster energy and the 
uncertainties involved in the treatment of both the tensor and the channel dependence of 
the realistic interaction forces. In this respect we explain how it is possible to believe that 
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the LOCV method is capable of giving good agreement with the results obtained by the 
more complicated schemes in which the many-body contributions have been taken in to the 
account. 



II. THE LOCV FORMALISM 

We consider a trial many-body wave function of the following form to evaluate the 
Rayleigh-Ritz upper bound to the ground-state energy: 

i^v = F<^> (1) 

where $ is a slater determinant of the plane waves of A independent nucleons, 

^ = Allexp{iki ■ fi) (2) 

i 

F is a A-body correlation operators which will be replaced by a Jastrow type i.e. 

F{l...A)=Sl[f{tj) (3) 

i>j 

and A and S are an anti-symmetrizing and a symmetrizing operators respectively. The 
cluster expansion energy is written as [13]: 



Em = T I =E, + E, + Es + --->Eo (4) 





H 




A < V 


'v \ 





where Eq is the through ground-state energy. In the lowest order, we truncate above series 
after £^2- 

The one body term Ei is just the familiar Fermi gas kinetic energy i.e. 

= (5) 

5 2m ^ ^ 

where kp — (|7r^p)^/^ and p, are the Fermi momentum and the density of uniform nuclear 
matter, respectively. The two-body energy E2 is: 

^2 = ^E<^i|V(12)|u>„ (6) 



and the "effective interaction operator" V(12) is given by tfie following equation: 

V(12) = -^[^(12), [VL, F{12)]] + F(12)F(12)F(12) (7) 

< ij \ O \ ij >a represent the antisymmetrized two-body matrix element taken with respect 
to the single-particle plane waves. 

The two-body correlation operator F{12) is defined as follows: 

m) = E f!f^Oj^\ij) (8) 

a,p=l 

a — {J, L, S,T} and operators 0^\ij) are written as: 

or-' = 1, + Isi,), (i - ^Sl,), SU (9) 

where 8(2 is the familiar tensor operator and is its analogous for the mixed A^-A channel 
(in case of the A-Reid potential) . The values of p is set to unity for L = and the spin- 
triplet channels with L ^ J ^ 1. But for L = J ± 1 it takes values of 2 and 3. Finally 
we have L — channels which couple the channel to the ^Dq channel (for the A-Reid 
potential) where we set p = 1 and 4. As in our previous works we also require that the 
correlation functions fa\fa^ fa^ ifa^) heal to the pauli function /p(r) (zero), 

fp{r) = [l-[h{kFr)]r'^' (10) 

where {Jj{x) are the spherical bessel functions) 

iW = ^ (11) 

The two-body nucleon-nucleon interaction V{12) has the following form, 

V(12)^J2V^^\ru)0^j (12) 

and they are taken from references [2,3]. 

The normalization constraint, < il^y \ i/^y >— 1, that we impose on the channel two-body 
correlation functions f^^ as well as the coupled and uncoupled differential equations, coming 
from Euler-Lagrange equations are described in references [1,2]. 



III. THE THREE-BODY CLUSTER ENERGY 



A. The state-averaged calculation 

In general the three-body cluster term expression in the energy expectation value has 
the following form [13], 

Es = + ) (13) 

where 



1 

ijk 



Ei'^ = -JJ:h^kW,, (14) 



Ei'^-4^E^^Jk (15) 



ijk 



and hik , Wij and Wijk are defined as, 



Kk =< I /i(12) I ik >a 

Wij =< ij I V(12) I ij >a (16) 
W,jk =< ijk I V(123) I ijk >a 

The two-body effective operator V(12) is given in equation (7) and {t{i) — — ^Vj) 

h{12) = F^(12)F(12) - 1 (17) 
V(123) = ^4(123) [t(l) + i(2)+t(3) , F3(123)] +adj. 

+ (Fj(123)y(12)F3(123) - V(12) + same for pairs 23 and 13) (18) 

Now by imposing the Jastrow approximation and ignoring the state and operator dependence 
of correlation functions, one can write the above three-body cluster energy E^ expressions 
in the following form : 
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where 

= E < I ^(^13) I >a< ij I V(ri2) I ij >a (20) 

ijk 

= E < I h{m)V{n2) I ijk >a (21) 

life 

Eil = ^ E < I ^(n3)V(ri2)/i(r23) I ijk >a (22) 
Tii = ^ E < I |^/'(n2) V2/i(ri2) • V2h{r2s) \ ijk >„ (23) 

ijk 

and the hole function h{rij) is defined as, 

Knj) = /^(r,,) - 1 (24) 

< ijk I O I ijk >a represent the antisymmetrized three-body matrix element taken with 
respect to the single-particle plane waves. 

Now, we define a state- averaged two-body correlation function [10] as 

.2. . ^ Ea,(2r+ 1)(2J+ 1)|[1 - {-l)'^-'^^'m^\x)af{x) 

E.,(2r + i)(2j + i)i[i-(-i)^+^+>r(x) 

and a state-averaged two-body effective interaction as 

^.^2) = E.,(2r+ 1)(2J+ l)i[l - (-l)^+^+^]V(')(12)ar (x) .^G) 
E.,(2T + 1)(2J+ l)i[l - (-l)^+^+^]ar(a;) 

where the ^^(a;), etc. are 

a(^^\x) = x\2J+l)-'[{J + l)Ij_,{x) + JIj+,{x)] (27) 

a(f)^(x) = x\2J+l)-'[JIj_,{x) + (J+ l)/,+i(x)] 
ai4)'(a;)=a;2j^(x) 

with 



Il{x) = 48 /' dzz\l - + ]:^')Jl{zx) (28) 
Jo II 

and 

m) = Hij) - 1 (29) 

Finally after doing some algebra, we get the following three-body cluster terms in terms 
h{ij) and V{ij) as, 
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{4 + h{kFr23)l(kFru)l(kFns) - (Pikpru) + f(kFru))} (30) 



{4 + hikFru)l{kFr23)l{kFns) - l\kFn2) - ll^kpr^^)} (31) 



Efi = ^ / dh.dh^dh, h(n,) h{T2,) V(ri2) 



{4 + -^l{kFn2)l{kFr2^)l{kFn^) - /'(fc^ris) - 2l\kFr2^)} (32) 

^(3) _ ^ [ ^3^^^3^_^3^^f^ 

8A 
2. 



Til = ^J d'nd'r2d'r,{—f\n,)V2Hn2) ■ V^Mr^s)) 

{4 + h{kFn2)l{kFr23)l{kFrn) - l\kFru) - 2f{kFr23)} (33) 



Then we can simply add the contribution of the three-body cluster energy, E3, to the two- 
body energy E2. 

— — f2) 

In the figure 1 and 2 we plot the contributions of various quantities such as E^, E^f^ , 
Efh^ Ezh = e'^h + Ef^, Efjlf^ , fg^^J^ and Efj^^ + Tg^^)^ against density in nuclear matter for 
the Reid and the A-Reid potentials respectively. As we expect, it is seen that the main 
contribution comes from E^h- Since the other terms namely E^hh and T^hh arc smaller than 
this term, i.e E^h with order of h{ij). It is also seen that the density dependence of E^ is 
mainly due to E^^. On the other hand the size of E^uh ^-nd T^uh ^-re roughly the same but in 
opposite sign beyond the nuclear matter density. So as we mention before, we can conclude 
that only E^h is important in the three-body cluster energy. 
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B. The state-dependent calculation 



Regarding our discussion in the section III-A , the main contribution of three-body cluster 

(2) (3) 

energy should come from E^,]^ and E^j^ . So in this section we explain how it is possible to 
write Ef^ and Ef^l in terms of f^\r) and V^\r) and calculate them exactly i.e. without 
making state-averaged approximation. 

By using the addition of angular momentum theorem [14], the partial wave expansion 
and various orthogonality relations, the explicit form of E^fl can be written as following, 

r 1. U 1. ^ T Q TT' T f Tfd/T'f 



{kikjkk}LSJT UJ'S'T' 



,(4vr) 
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E (34) 

MjMj,\Ml Ms Mj) \Ml' Mg' MjJ 

where the large parentheses are the familiar Clebsch-Gordan coefficients and 



h^s-^T{x) =< LSJT I h{x) I LSJT > V^^^^iy) =< LSJT \ V{y) \ LSJT > 

kik — 2^^* -^ik — ki -\- k). Kij = ki -\- kj 

cr's and r's stand for particles spin and isospin respectively. The following identity, 

E M^jh iPx)jh iQy)jl2 iPy)Yhrm iP)Yl2m2 iP)Yllm2 (^) = 

hl2mim2 



(47r) 
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^ J dVL^dVLy exp («(f - y) ■ {q- p)) (35) 



is used to check the vahdity of equation (34). To see this matter clearly we assume that 
V and h are state- independent, then, it is possible to perform all of the summations and 
we show that one can reach to our previous equations for the state-independent three-body 
cluster energy (equation (30-33)). 

In figures 3 and 4 we plot by using the state-dependent correlation functions and 
effective potentials i.e. fa\'^) ^a'*(^) for Reid and A-Reid Potentials respectively. 
Ef^{C) shows the same calculations but only with the central parts of f^\r). Effl{NC) 
gives the estimate of non-central parts i.e. El^-Ely'{C). Prom these figures, as one should 
expect, it can be concluded that the main contribution in E^y^ comes from the central parts 
of correlation functions and with very good approximation we can ignore the non-central 
parts of energy contributions in the three-body expressions (this is not a good approximation 
in the two-body level). 

The exact calculation of E^^l is very comphcated, but according to the above argument 
we can take into the account only the central part of correlation functions and make it 
possible to have an estimate for this term as well i.e. 

1 



Etac) = -]4 E < I /i(13)V(12) I ijk >a (36) 



ijk 

We start by inserting a unit operator J2imn I ^'^^ >< ^'^^ I between the two-body operators 



f{ij) and V{ij) in above equation. 

£;^J(c) = l <ijk\h{13)\lmn><lmn\V{12)\ijk>a (37) 

ijk,lmn 

Now by expanding the anti-symmetrized ket we get 9 terms which are very similar when we 
want to calculate them. We explain how it is possible to calculate one of these terms i.e. 

^ < ijk I h{13) I Imn >< Imn \ V(12) | ikj — kij >= 

ijklmn 

Y.<ik\ /i(13) I Ij >< Ij I V(12) \ik-ki>^ ^ EE/ / dxdy 



ijkl {kikjkkki} LSJT L'J 

03 



(^?/)'m(^^^^^)((-l)^+^+^ - l)h^''^{x)V^''''^{y) (38) 



JL{kikx) Jl ikjx) Jl' {kiky) Jv {hjV) 
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( L S J V ( L' S J' Y ^ ^ ^ _ 
MjM,jMs \Ml Ms Mj ) \ Ml, Mg Mj, J 

In figures 5 and 6 we plot the various quantities such as E^i^{C), E^^iC) and 
E^h{C)^Effl{C) +Effl{C) for the Reid and the A-Reid potentials respectively. ^S^C), 
E^f^(C) and E^h{C) are the same terms but with the state-averaged approximation of sec- 
tion III-A. 

Finally in figures 7 and 8 we plot the total contribution of E^h{C) by presenting the 
exact and the state-averaged calculation for the Reid and the A-Reid interactions. The size 
of £"3/1 (C) is less than state- averaged E^hiC) and this difference becomes sizable beyond 
the nuclear matter saturation density. This shows that at large densities the state-averaged 
approximation is no longer valid. Our final results for the binding energy of nuclear matter 
per nucleon are given in table 2 by adding the three-body cluster energies to our LOCV 
calculations. 

IV. DISCUSSION 

We have developed a method to make a very good estimate of the three-body cluster 
energy by using the state-dependent hole functions and effective interactions. We found that 
the size of the three-body cluster energy is much more smaller than our previous estimate 
using the state-averaged approximation. The smallness of the three-body cluster energy 
indicates that we are very close to the exact results at LOCV level and our cluster expansion 
series converges very rapidly. Our result with the A-Reid interaction is comparable with 
those calculations in which the three-body interactions [5-9] have been taken into the account 
such that they could reproduce the empirical properties of nuclear matter [5-9] . Presumably 
in the lowest order limit the cluster expansion for the exact two-body radial distribution 
function, g{r)^ is automatically convergent and, provided that the lowest order energy has 
been efficiently minimized, the replacement of g{r) by its lowest order approximation, 
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'^'\r) = {1 - 



(39) 



is then a very good approximation. In this case the integral constraint (normahzation 
constraint) , 



plays a crucial role at all densities in forcing the correlation function to heal and in restricting 
the size of the wound in the correlated two-body wave function. 

It is interesting to note that it may also be the inclusion of the tensor and spin-orbit 
forces and the incorporation of tensor correlation which allows us to deal accurately with 
the contributions to the energy of nuclear matter from all of the states of reahstic potentials 
and to obtain such good agreement with far more ambitions and complicated calculations 
such as VHC. 

Whether or not the above arguments are quantitatively correct, we still have to face that 
for p < 0.3/m~^ our nuclear matter results agree extremely well with those of VHC, which 
are hopefully reliable upper bounds on the energy up to this density and differ considerably 
from the accepted experimental equilibrium values. So we can conclude that only by taking 
into accounts the internal degrees of freedom of nucleons (such as A(1234)) [2], the three- 
body forces [5-9] or the relativistic effect [16] it may be possible to remove this discrepancy 
between the theoretical calculation and the empirical prediction. 

In this context we remark again that the LOCV method because of (1) its agreement with 
VHC which includes the many-body cluster contributions, (2) the smallness of three-body 
cluster energy and (3) its great simplicity, is a useful tool in the study of other properties 
of nuclear matter and finite nuclei. 




(40) 
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Figure captions : 

Figure 1: Various terms of the three-body cluster energy (MeV) versus density {fm~^) with 
the state-averaged correlation function and the effective interaction in nuclear matter ac- 
cording to the equations of section III-A for the Reid interaction. 

Figure 2: Same as figure 1 but for the A-Reid interaction. 

Figure 3: The plot of (MeV) versus density (/m ^) by doing an exact calculations 
with the Reid potential. £;^J(C) and Efjl{NC) are the central and non-central parts of this 
term, respectively. 

Figure 4: Same as figure 3 but for the A-Reid potential. 

Figure 5: Comparisons of E^h^C) (MeV), etc. with E^hiC) etc. versus density {fm~^) 
for the Reid potential 

Figure 6: Same as figure 5 but for the A-Reid potential. 

Figure 7: The plots of the state-averaged E^h and its exact central contribution, E^h{C) 
(MeV) versus density {fm~^) for the Reid interaction. E^h^C) is the central part of state- 
averaged contribution 

Figure 8: Same as figure 7 but for the A-Reid potential. 
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Table 1. The saturation energy, the saturation density and the incompressibihty of nu- 
clear matter for different potentials and the many-body methods (see different abbreviation 
and references on page 2). 



Potential Method Author po{fm-^) E(MeV) /C(MeV) 
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340 




LOCV 
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238 
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Table 2. The saturation energy and the saturation density of nuclear matter by adding 
the three-body cluster energy according to the text. 



Reid 


A-Reid 


Eo E-, + E,, Ei + Ey, 


E2 E-, + E, Eo + Es 


Poifm-^) 0.294 0.216 0.295 
E(MeV) -22.83 -22.24 -23.362 


0.258 0.208 0.257 
-16.74 -15.597 -16.721 
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